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If dark matter in tiie galactic halo is composed of bosons that form a Bose-Einstein condensate then 
it is likely that the rotation of the halo will lead to the nucleation of vortices. After a review of the 
Gross-Pitaevskii equation, the Thomas-Fermi approximation and vortices in general, we consider 
vortices in detail. We find strong bounds for the boson mass, interaction strength, the shape and 
quantity of vortices in the halo, the critical rotational velocity for the nucleation of vortices and, in 
the Thomas-Fermi regime, an exact solution for the mass density of a single, axisymmetric vortex. 
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I. INTRODUCTION 

The observational evidence for the existence of non- 
baryonic dark matter is impressive. Precision measure- 
ments of anisotropies in the cosmic microwave back- 
ground reveal a universe composed of 4% baryons, 22% 
non-baryonic dark matter and 74% dark energy [1]. Ob- 
servation of galaxies and galaxy clusters from rotation 
curves [2], gravitational lensing [3] and X-ray spectra [1] 
has lead to a picture in which galaxies are composed of a 
luminious galactic disk surrounded by a spherical galac- 
tic halo of dark matter which comprises roughly 95% of 
the total mass of the galaxy [5] . 

Various explanations for dark matter have been pro- 
posed. Since we are interested in a particle explana- 
tion we will mention the two most prominent proposals 
from particle physics: WIMPs (weakly interacting mas- 
sive particles) [6] and axions [7]. WIMPs have a mass 
around a TeV and interact via the weak force. Weak in- 
teractions automatically lead to the correct relic density 
[H]. A well known WIMP candidate comes from super- 
symmetry and is known as the lightest supersymmetric 
partner (LSP). In many scenarios this is the lightest neu- 
tralino, which is a linear combination of the fermionic su- 
perpartners of the Higgs scalars and neutral electroweak 
gauge bosons. Axions are scalar bosons first proposed 
to solve the strong CP problem in QCD and may be 
produced, for example, by an initial misalignment of its 
vacuum angle in the early universe j5]. 

An intriguing possibility is that if dark matter is com- 
posed of bosons (which includes, but is not limited to, 
the axion |18j ) then at sufficiently low temperatures 
the bosons will form a Bose-Einstein condensate (BEC). 
Galactic halos could then be gigantic BECs. Such an 
idea, that dark matter may be in the form of a BEC, 
has been proposed by a number of authors [5HT1] and 
has been shown to alleviate known problems with the 
standard cold dark matter picture, such as cuspy cen- 
tral densities and an abundance of small scale structure 
in halos TT. The standard approach to treating a non- 
relativistic BEC is to use the Gross-Pitaevskii equation, 
also known as a nonlinear Schrondinger equation, which 
includes a self-interaction term and a confining, trap po- 



tential. For gravitating BECs the trap potential is simply 
the Newtonian gravitational potential, so that the Gross- 
Pitaevskii equation is coupled with the Poisson equation 

[13 [Hill]. 

A rotating BEC is known to give rise to a lattice of 
quantized vortices [551 US] ; without which the condensate 
would be irrotational. The formation of vortices leads to 
lower energy states and thus is energetically preferred. 
If we are imagining the galactic halo to be made up of 
bosonic dark matter in the form of a BEC the rotation 
of the halo may give rise to the formation of vortices 

(laiisHni. 



If vortices form in the halo, they may affect the mass 
distribution and other properties. In this paper we con- 
sider vortices in BEC dark matter in detail. In the 
next section we review the Gross-Pitaevskii equation, the 
Thomas-Fermi approximation and vortices. These equa- 
tions will be the framework in which all of our work is 
done. In section IIIII we consider the condensate as a 
whole, populated by a lattice of vortices, and make a 
model independent argument that for a stable conden- 
sate with vortices the boson should have a mass centered 
around m ~ IQ-^" - IQ-^s kg ~ IQ-^* - 10-^2 eV. In 



section IV we analyze a single, axisymmetric vortex and, 
in the Thomas-Fermi regime, find an exact solution for 
the mass density. Our results lead to a comprehensive 
picture for both the condensate and vortices, which we 
present in sectionjvj Using our vortex solution, in section 
VI we compute the critical rotational velocity of the halo 
such that vortices would form and find that it is likely to 
have happened. We conclude in section |VlT| 



Although much of this paper is theoretical, it is im- 
portant that we connect our results to data. When do- 
ing so we will use the following numbers from M31 (the 
Andromeda galaxy) [191 , which we take to be not unrea- 
sonably characteristic of dark matter halos: A total mass 
of M ^ IQ^'^Mq ~ lO''^ kg, an average radius of roughly 
R ~ 100 kpc ^ 10^^ m and an approximate rotational 
velocity of 125 km/s at 30 kpc. From this we estimate 
an average mass density of pavg ^ kg/m"^ and a ro- 

tational angular velocity for the halo of '-^ 10~^^ rad/s. 



II. GRAVITATING BOSE-EINSTEIN 
CONDENSATES AND VORTICES 

In this section we give a general introduction to the 
Gross-Pitaevskii equation, the Thomas-Fermi approxi- 
mation and vortices as they apply to gravitating BECs. 
We give considerable attention to these equations as all 
of our results are derived from them. 



A. The Gross-Pitaevskii Equation 

We assume that dark matter is composed of scalar 
bosons, of mass m, that have undergone a phase tran- 
sition towards Bose-Einstein condensation. To describe 
the BEG we use the standard symmetry-breaking mean- 
field approach, which is expected to be valid for systems 
with a sufficiently large number of particles and at tem- 
peratures far bellow the BEG transition temperature. 
In this approach one associates an order parameter, or 
macroscopic wave function, ipi^r), with the bosonic field 
and ignores high-order correlations due to bosonic quan- 
tum field fiuctuations. The BEG can then be described 
by the time-independent Gross-Pitaevskii energy func- 
tional: 



1, 



1 



The first term represents the kinetic energy and the sec- 
ond term coUisional energy due to repulsive contact in- 
teractions with strength 



ATTlT?a 



(11.2) 



where a (> 0) is the s-wave scattering length. The final 
term is the gravitational potential energy, with gravita- 
tional potential 



Vair) = -Gm J 



which satisfies Poisson's equation: 



(11.4) 



The condensate wave function, ^, obeys the Gross- 
Pitaevskii equation 



V' + Vo\M' + mVG]^ = t^^, (II.5) 



2 m 



which is obtained either by variation of E[ip] with re- 
spect to ip* under the constraint that the total number 
of particles, N, is conserved (i.e. J d'^r|-!/)p — N), or 
equivalently by variation of the grand-canonical thermo- 
dynamic potential, Eltp] — fiN, with respect to where 
/i is a Lagrange multiplier representing the chemical po- 
tential. Our system is then described by the coupled 
differential equations (II.5) and (II.4|. 



It is very common, and useful, to write the wave func- 
tion in terms of its modulus and phase: 



V;(r) = |7/^(r)|e*^W, 



(11.6) 



both of which are real, describing a condensate with local 
velocity 

vir) = -VS{r) (II.7) 
m 

and mass density p{r) = m\Tp{r)\'^ with normalization 

d^rp = M, (II.8) 



where M is the total mass of the condensate. 

In terms of these hyd ro dynamic variables the Gross- 
Pitaevskii equation (II.5) and the Poisson equation (II.4) 
become 



= V • (pv) , 

V^Tp 1 Vo 
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2 



= V^Vg - AttGp, 



(11.9) 

p + mVc-M, (11.10) 
(11.11) 



and the energy functional (II.l) takes the form 



(II.l) E[p,v] 



^0 2 



Vg 



2to2' 



-p 



2 

(11.12) 

where (II. 9) and (11.10) can be obtained either from 



equating the imaginary and real parts of the Gross- 
Pitaevskii equation (II. 5 I or by minimizing (11.12) with 
respect to v and p. 



Aside from the first term in (11.10) (which is called the 
quantum pressure term), Eqs. (II. 9) and (II. 10 1 make 



up the (time-independent) hydrodynamic description of 
the condensate since (II. 9 1 is the continuity equation and 



(11.10) is the Euler equation from classical fluid dynam- 



(II-3) ics. The quantum pressure term in (11.10) can be traced 



to the first term in (11.12) which originates from the un- 



certainty principle and hence cannot find its analog in 



classical physics. In contrast to the other terms in (11.12) 



this quantum kinetic energy contains the gradient of the 
mass density p{r). As a result, as the number of par- 
ticles in the condensate, or equivalently the size of the 
wave function, increases the quantum kinetic energy be- 
comes negligible compared to other energy contributions 
except near boundaries of the condensate. Neglecting 
the quantum pressure term (i.e. the quantum kinetic en- 
ergy) is known as the Thomas-Fermi approximation. It 
is important to realize that in our case gravity does not 
constitute an external trap fixed, a priori, by an exter- 
nal agent. Instead it depends on the wave function and 
hence must be determined self-consistently with the help 
of the three coupled equations (II.9), (11.10) and (11.11 1. 



In analogy to the ordinary Thomas-Fermi regime which 
is defined to be dominated by the balance between the 



2 



external trap potential and the repulsive s-wave scatter- 
ing, we define, as in [20] . the gravitational Thomas- Fermi 
regime as the regime dominated by the balance between 
the gravitational field created by the condensate parti- 
cles and the repulsive s-wave scattering. In this regime 
the quantum pressure term may be ignored and (11.10) 
reduces to 



1 



= —rav H p + raVc — M 

2 m 



(11.13) 



from which we determine the coefficient in (11.161. Since 



we are considering the case of rotating around the z-axis 
we find 



5{r ~ ro)z. 



(11.21) 



Vortices with k ^ 1 are generally unstable and decay 
into fc = 1 vortices [23j. As such, we will confine our 
attention to fc = 1 vortices. 



and the energy functional (11.12) to 



E[pM= I 



(11.14) 



Vortices 



To motivate the need for vortices note from (II. 7 1 that 



since the velocity of the condensate is defined in terms of 
a gradient, its curl vanishes: 



V X i; = —V X (V5) = 0, 
m 



(11.15) 



suggesting that a BEC is irrotational. Since BECs are 
well known to rotate v must contain a singularity, so as to 
evade a vanishing curl, generalizing the above expression 
to 



V X D ^ S{r — rg). 



(11.16) 



This singularity is the vortex core, centered at rg, and 
represents the absence of condensate (i.e. the mass den- 
sity is zero inside the vortex core). 

We imagine the vortex to be cylindrical in shape and 
adopt cylindrical coordinates (r'j_, z). The wave func- 
tion takes the form 



^(r) = |7/^(r^,z)|e^'=^ 



(11.17) 



which we refer to as the vortex ansatz and which de- 
scribes a fc-quantized vortex. The 27r-periodicity of the 



wave function requires k to be an integer. From (II.7) we 
find 



hk 7 

V = 4>, 

■mr± 



(11.18) 



and thus this is a vortex rotating around the z-axis. 
The circulation is given by 



K ^ (h d£ ■ V = 



2Trnk 



(11.19) 



which may also be evaluated using using Stoke's theorem, 
2TThk 



J dA-{V xv), 



(11.20) 



III. BOUNDS ON THE BOSON MASS 



Both the energy functional (II.l) and the Gross- 



Pitaevskii equation (II.5) were written in a non-rotating 
frame. Since we are considering rotating galactic halos, 
it is natural to work in a frame that is rotating with the 
same angular velocity as the halo and where the poten- 
tial energy is time-independent. Such a transformation 
(to a frame rotating with constant angular velocity fl) 
modifies the energy functional by 



Em ^ Em 



(III.l) 



where L = rx {h/i)\J is the angular momentum operator. 



so that (11.14) becomes 



E[p,v] = J 



lv'-v{nxr) + ^p + VG 



(111.2) 

Without loss of generality we assume the frame to be 
rotating around the z-axis so that the constant angular 
velocity fl points in the positive z-direction. We adopt 
cylindrical coordinates (rj_,(f),z) and define 



ft X r = Vlr^ij). 



(III.3) 



We may rewrite (III. 2) as 



Em - 



I , .2 1 

-{v^vn) - -V 



Vn 



(ni.4) 

showing that the value of v which minimizes the energy is 
precisely v — vq. Now, V x vq = 2fl, while from (11.21 ) 
W X V = {2nhk / m)5{r — ro)^. Thus it is not possible 
for a condensate to acquire the solid-body rotation v = 
Vfi. Instead it will nucleate a lattice of uniformly spaced 
vortices in an attempt to mimic solid-body rotation |23) . 
In this sense v ^ vq. 

If the BEC has constant angular velocity andis pop- 
ulated by a lattice of A'vor vortices then from (11.21 ) 



V X V 



27rfi 
m ^-^ 

i=l 



S{r 



(III.5) 



where we have assumed all vortices have fc = 1 as men- 



tioned at the end of section II B The circulation (11.19) 
gives 

(III.6) 



KN = f M-V ^ 27rf7i?BEc = 2rjABEC, 



3 



where the path was taken around the entire condensate 
and thus i?BEC is the radius of the condensate and ^bec 



its area. Using Stokes theorem, as was done in (11.20), 
and (III. 5 1 we find also for the circulation 



J dA- {\/ xv) 



2TTh 



(III.7) 



Equating our two results we determine the areal density 
of vortices [211 : 



A 



EEC 



(III:, 



The inverse of this quantity is the area per vortex. As- 
suming it to be circular with radius i?vor, so that ttR^^^ — 



1/nvor, we find 



h 



(III.9) 



With this equation we can find a rough lower bound on 
the boson mass. The size of the vortex must be smaller 
than the condensate: i?vor < -Rbec- Since we are assum- 
ing that dark matter in the galactic halo is in the form 
of a BEC, then with the numbers from the end of the 
introduction wc find 



-^BEC^ 



10"^° kg ~ 10-^*eV. 



(III.IO) 



This is approximately the mass many authors expect 
the boson to have for BEC dark matter El [H [H 
1161 117) . The reason is rather simple and can be under- 
stand from a simple qualitative argument [TTJ [26l l30] : 
For an ultralight dark matter particle with a Compton 
wave length of astrophysical size, the uncertainty prin- 
cipal can suppress the formation of excess small scale 
structure seen in the standard cold dark matter scenario. 
The masses considered obey 



m < 10"^^ kg - 10"^^ eV 



22 



(III.ll) 



(despite the small mass, BEC dark matter behaves like 
cold dark matter [TTJ [551 HI])- In combination with 
(III.IO I we find that for a galactic halo of BEC dark 



matter with vortices there is surprisingly little param- 
eter space available for the boson mass. Our model inde- 
pendent arguments have lead to a boson mass centered 
around 

m - 10"^" - 10"^^ kg - 10-2"* - 10^22 eV. (III.12) 

IV. A SINGLE AXISYMMETRIC VORTEX 

In this section we solve for the mass density of a single 
vortex centered on the axis of rotation. The condensate, 
in the Thomas-Fermi regime, is described by the Gross- 
Pitaevskii equation (11.13), 



1 2 Vo 

= -mv H p + mVc — 

2 TO 



(IV.l) 



and the Poisson equation (11.11), 

V^Vg = 47rGp. 



(IV.2) 



We would like to find a vortex solution to these coupled 
differential equations. To do so will use a combination of 
spherical coordinates (r, 0, 6) and cylindrical coordinates 
(rj^,0, z) and assume the solution to be of the form of 
the vortex ansatz (11.17) with fc = 1: 



(IV.3) 



where p{r,0) — p(r±,z) = m\ip\'^ is the mass density. 
From (II.7) the velocity of the condensate is 



V = 



-0 



and the Gross-Pitaevskii equation becomes 



= 



2mr\ 



H p + mVG — P- 



(IV.4) 



(IV.5) 



Taking the Laplacian of this equation and combining it 
with the Poisson equation (IV.2) we arrive at 



= 



Vn 



-P 



Vo ri' 



(IV.6) 



Our system is described by this single equation. 

It is convenient to write this equation in terms of di- 
mensionless variables, O and £ — {f, f±,z}, defined by 



\ Vn 



Gm^ 



(IV.6 1 may now be written 

= v^e 



(IV.7) 
(IV.8) 

(IV.9) 



It is important to stress that while the original Gross- 
Pitaevskii equation is nonlinear, the equation for the 
mass density under the Thomas-Fermi approximation is 
linear. This allows us to write the solution to (IV.9), 



(IV.IO) 



as a superposition between 0/j (f , 6) — the solution to 
the homogeneous part of the differential equation — and 
0p(^_l) — the particular solution originating from the 
1/f^ term. The homogeneous part is the Helmholtz 
equation with solution 



E 

1=0 



aiji{f) + bm{f)] Pi {cos d), (IV.ll) 



4 



where ai and bi are arbitrary constants, ji and yi are 
spherical Bessel functions and Pi are Legendre polyno- 
mials. The particular solution is 



1 



1/2,3 
1,1,2,1/2 



1,2,1 



(IV. 12) 



where Jq and Yq are Bessel functions of the first and 
second kind and 



/"im.n 



ai, 



I Op 
, bn 



(IV.13) 



is the Meijer G- function [57], and is well approximated 

by 



1 



(IV. 14) 



In order to fix the constants in (IV. 11 1 we note that 
since the influence of the rotational energy (the 1/f^ 
term in (IV.9|) is vanishingly small in the large f± limit. 



the vortex solution must match the vortex free solution 
for large r± in the Thomas- Fermi limit. This observa- 
tion translates our task to finding all constants for the 
vortex free solution Qh{r, 0). The vortex free solution is 
well-known in the literature }28| , having the form 



(IV.15) 



which is obtained from (IV.ll) by requiring it to be 



spherically symmetric and finite at the origin which sets 
all constants to zero except ag = Go- 
The vortex solution is then 



1 



(IV.16) 



and the number of vortices in the halo is roughly 
RlaJR^ ~ 1 - 100, where R - lO^i m is the radius 
of the halo from the end of the introduction. 

To analyze the shape of the vortex we use the solution 
(IV.16). We expect it be roughly spherical with a thin, 
cylindrical vortex core along the z-axis. Since the vortex 
core is along the z-axis, to determine its size we analyze 



our solution ( IV.16| for small fj^. The region of negative 
mass density is taken to represent the core itself, the 
mass density being negative because the Thomas-Fermi 
approximation is invalid in this region. We will use the 
standard technique of cutting off the solution, removing 
the region of negative mass density and setting this region 
equal to zero. Thus, the vortex core is defined by the 
smallest values of f j^c which satisfy 



On 



(V.3) 



The z-dependence, around z = and for small r_Lc, is 



1 



2ve 



The radius is smallest in the z = plane: 



1 



(for z = 0). 



(V.4) 



(V.5) 



In the unsealed unit system this equation simply corre- 
sponds to the healing length, R^ = \J m/Sirapa, where po 
is the mass density at the center of a vortex free conden- 
sate. 

The edge of the vortex is taken to be the first point, 
outside the core, where the mass density is once again 
zero [131 im HH] ■ Away from the upper and lower caps 
(in the z-direction) the — l/4f^ term is negligible and the 
edge of the vortex is defined by 



eo 



smrv 



0, 



(V.6) 



V. CONDENSATE AND VORTEX STRUCTURE 

In this section we compile and analyze our results in an 
attempt to say as much as possible about the structure 
of the condensate and the vortices which reside in it. We 
will form a comprehensive and detailed picture. While 
we are using the numbers from the end of the introduc- 
tion, and thus our results are specific to the Andromeda 
galaxy, we expect them to hold more generally. 

In section Hill we found that the boson mass is centered 
around 



10"^° - 10 



-58 



kg ' 



10^2^ - 10"^2 eV. 



(V.l) 



Returning to (|III.9|) and using these masses, the size of 

10^0 



the vortex is around 



10^^ m 



(V.2) 



with solution fvor = tt. Thus, the vortex is roughly spher- 
ical in shape with radius tt or, in the unsealed unit sys- 
tem. 



(V.7) 



In figure [T] we plotted both the edge of the vortex and 
the edge of the vortex core for Go = 1. The figure on the 
top shows the spherical shape of the vortex. The opening 
to the vortex core is visible. The figure on the bottom is a 
cross section showing the core lying along the z-axis and 
cylindrical in shape. The core can be made arbitrarily 
thin by increasing Qq. The figure is in good agreement 
with the actual vortex state only in the Thomas-Fermi 
limit where the healing length is much less than the size 
of the condensate, or equivalently i?c/i?vor — ^'^c/t^ — 
l/27rye;;< 1. 
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FIG. 1. The edge of the vortex and the edge of the vortex 
core are plotted for Go = 1. The figure on the top shows 
the spherical shape of the vortex with the opening to the 
vortex core visible. The figure on the bottom is a cross section 
showing the cylindrical shape of the core. The core can be 
made arbitrarily thin by increasing Oo- 



We may use (V.7) to determine the chara cteristic size 
of the scattering length a. Comparing (V.7) to the char- 
acteristic size of vortices given by (III. 9) we find that the 
bosons should obey 



G 
Ml 



(V.8) 



In lieu of our expected boson mass ( V. 1 ) the scattering 



length is expected to be centered around 



10-80 _ 10 



-76 



(V.9) 



Determining the characteristic radius of the vortex core 
from (V.5) requires determining Gq. To determine we 



II 




5f 0.001 



10-" 



FIG. 2. The radius of the vortex core, in the z = Q plane, 
as a fraction of the radius of the vortex is plotted versus the 
boson mass. The various lines correspond to different scatting 



lengths. From upper left to lower right: 



10" 



10" 



10" 



m, 10" m, 10"'^ m and 10"'^ m. 



might try fitting our solution to data, but since we expect 
the galactic halo to be populated by many vortices it is 
not simple to fit the solution for a single vortex to data 
for an entire galactic halo. However, we may determine 
00 if we make the following assumption: We assume that 
the average mass density of our vortex solution is char- 
acteristic of the average mass density of the galactic halo 
in which it resides. In the Thomas-Fermi regime the core 
has a negligible affect on the total mass and we find 



Pavg 



(4/3)7ri?3^, 7o 



dr p{r) = 



2TT^h^a 



From (V.5) the vortex core size is then 



(. 



1/2 



(V.IO) 



(V.ll) 



In figure [2] we have used (V.ll) and, from the introduc- 
tion, pavg ~ IQ-^'^ kg/m^ to plot the radius of the core 
in the z = plane as a fraction of the radius of the vor- 
tex versus boson mass. The different lines in the figure 
correspond to different values of the scattering length, a, 
as explained in the caption. 

Our results have lead to a comprehensive picture of 
the condensate and vortices. To reiterate, if dark matter 
in the galactic halo is composed of bosons in the form 
of a BEC that has formed vortices then we should ex- 
pect the boson mass to be centered around m ^ IQ-^^ — 
10-58 kg IQ-^"^- 10-22 eV, the scatting length around 



10 



-80 



— 10 m, the radius of the vortex around 



1920 — 102^ m, which corresponds to roughly 1 — 100 vor- 
tices in the halo, and a vortex core radius as displayed in 
figure [2j The vortex is expected to have a shape that is 
nearly spherical with a thin, cylindrical vortex core along 
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the axis of rotation, as displayed in figure [T] and a mass 
density described by (IV. 16). 



VI. 



CRITICAL VELOCITY FOR VORTEX 
FORMATION 



and can be conveniently calculated entirely in the labo- 
ratory (unprimed) frame. 

In the following subsections we calculate the necessary 
energies and approximate the critical angular velocity at 
which vortices form. 



The vortex solution in section IIVI was obtained in a 
non-rotating frame. The fact that the galactic halo ro- 
tates with angular velocity n motivates the question of 
whether a given value of O can sustain a single vortex 
or an array of vortices. In this section we calculate the 
critical angular velocity of the condensate at which vor- 
tices begin to form. To understand how [29] recall from 
(III.ll that when moving to a frame rotating with con- 
stant angular velocity J7 the energy transforms as 



E' = E- 



(VI.1) 



In the rotating (primed) frame we calculate the energy 
for a condensate with a single (fc = 1) vortex, which we 
call E[(il.). We also calculate the energy for a vortex 
free condensate (i.e. the fc = configuration), which we 
call E'q. In the latter case, since there is no rotation, 
the energy is independent of and may be calculated in 
either frame {E'q = Eq). The difference in energy. 



AE'{n) = E[{n) -e'q = E[{n) - Eq 



(VI.2) 



is the energy required to form a vortex at angular ve- 
locity i7. The formation of vortices will occur when it 
is energetically preferred, i.e. when /S.E'{Vi) < 0. The 
critical velocity for vortex formation is then defined as 



AE\n,) - 0, 



(VI.3) 



and vortex formation will occur for > flc- 

For an axisymmetric trap potential, such as the gravi- 
tational potential, there is a particularly nice formula. In 
this case, if we take the condensate to be rotating around 
the 2;-axis the wave function is an eigenfunction of the 
operator and an axisymmetric vortex has total angular 
momentum 



Mh 
m ' 



(VI.4) 



where M is the total mass of the condensate. Thus, from 



( VI.l I we have 



E' = E -n 



Mh 



m 



and 



AE'{n) = E[{n) ~ E^{n) = Ei-Eo-n 



Mh 



The critical angular velocity is then given by 

Ttl 

rip — (E-\ — En] 

Mh^ 



(VI.5) 



(VI.6) 



(VI. 7) 



A. Vortex Free Energy 

A vortex free condensate (i.e. the k = configuration), 
has, from (11.17) and (II. 7), ip{r) = \il'{r±,z)\ and no 
rotation, v = 0. To calculate the energy we use the 



energy functional (|II.14|): 



Vo 



The corresponding Gross-Pitaevskii equation is 
Vo 



-p + mVa- 



(VI.8) 



(VI.9) 



While this mass density, p, differs from the mass density 
when there is a vortex it must satisfy the same normal- 
ization: 



(fr M, 



(VI.IO) 



where M is the mass of the condensate, which doesn't 
change as a vortex is formed. 

Multiplying (VI.9) by p and integrating we find 



2m 
~M 



-En 



(VI. 11) 



We assume that p, is constant and thus the right hand side 
of (VI.9 1 must be independent of position. This means we 



are at liberty to evaluate the right hand side at any point 
inside the condensate we'd like. We choose to evaluate 
it in the z = plane at the edge of the condensate, 
rj_ =^ -RbeC: where the mass density is zero. Combining 
the two equations for p we have 



Eo = ImVg{Rbec)- 



(VI.12) 



We will comment on this term below after looking at the 
contribution from the vortex. 



B. Single Vortex Energy 

For a single, axisymmetric vortex we use the vortex 
ansatz (11.17), with A: = 1, which has velocity s quared 
v'^ = hy m r ^. The energy is still obtained from (|II.14[), 



El 



3 / 1 2 , ^0 1 

r I —V H 

2 2m2 2 



(VI.13) 



and the Gross-Pitaevskii equation is 
1 

2 TO 



p = -TOV^ -t- — p -h mVa- (VI.14) 
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While this mass density, p, differs from the mass density 
for the vortex free condensate considered above, it must 



satisfy the same normahzation (VI. 10) with the same 
mass, M . 

Multiplying (VI. 14) by p and integrating we find 



m 
M 



2Ei 



1 



(VI.15) 



As before the right hand side of (VI. 14) must be inde- 



pendent of position and we choose to evaluate it in the 
z = plane at the edge of the vortex, r±^ = i?vor, where 
the mass density is zero. Combining the two equations 
for p we have 



1 



1 



El = -MVg{R.o.) + -Mv'{R^,,) + Ei^„, (VI.16) 



where 



(Pr 



(VI.17) 



Ei T, appears to be impossible to evaluate analytically. 
Building upon the intuition gained from figure [T] we es- 
timate the integral by approximating the vortex solu- 
tion (IV. 16 1 with the vortex free solution ( IV.15[ ) and an 
empty cylindrical core along the z-axis with constant ra- 
dius rj^c- This estimation gives an upper bound on the 
integral. Under this approximation Ei ^, amounts to, in 
spherical coordinates. 



Gm f 
dfdOdcj) sin0- 



h 

8^ 



2-77 



TT— tan ^(r±c/'^) 



tan^^(r^c/^) ^ic/ ^ 



dr ■ 



(VI.18) 



To leading order in rj_c this evaluates to 
Mh 



El, 



Si(7r)ln ( 2— | +r 



(VI.19) 



where Si(7r) « 1.852 is the sine integral evaluated at tt, 

• w -2.658, (VI.20) 



32 -3! 52 • 5! 



and 



c,0 



1 

8M 



GmS 1 



4 mi?|EC ' 



(VI.21) 



where to obtain the second equality we used ( V.5 ), ( V.7 1, 
( V.ll I and pavg = M/ {4:/3)ttRI^^ and replaced i?vor with 



i?BEC since we are assuming that a single vortex is form- 
ing over the entire condensate. 
Putting everything together: 



C. Critical Angular Velocity 

The critical angular velocity for the formation of vor- 
tices is 



Si(7r) In 2 



where from (V.ll I 



1 



''_Lc = 



2 \2TT^h?a^p^ 



1/2 



- + r , 

T 

(VI.23) 
(VI.24) 



and Qc,o is given in (VI.21 1 and T in (VI.20). In figure 
|3]we plotted il,c/^c,o versus the radius of the core in the 
z = plane as a fraction of the radius of the condensate. 
One can see that ric/^lc,o £ 0(10). With the numbers 
from the end of the introduction and our range of boson 
masses (V.l) we have il.c,o ~ 10~^^ — 10~^^ rad/s and 
thus there exists large portions of parameter space where 
the Andromeda galaxy would have formed vortices. 



El ^-MVaiR^or) 



1 Mh, 



Mh, 



TT m 



m 



Si(7r) In 2 



(VI.22) 



It is not difficult to show that the first term, like its coun- 



terpart in (VI. 12), is roughly hy^ Gm/a^, without any 
factors of f±c m the denominator, and thus negligible. 
The second line above is then the dominant contribution. 



VII. CONCLUSION 

We considered the possibility that dark matter in a 
galactic halo is composed of bosons in the form of a Bose- 
Einstein condensate. Rotation of the halo is likely to nu- 
cleate vortices and we studied such vortices in detail. We 
found, in the Thomas-Fermi regime, an exact solution for 
the mass density of a single, axisymmetric vortex and, by 
fitting our equations to astrophysical data from the An- 
dromeda galaxy, developed a comprehensive picture for 
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10"'' 0.001 0.01 0.1 1 



rj.clrA.mc - RJRbec 

FIG. 3. fic/fic,o is plotted versus the radius of the core in the 
= plane as a fraction of the radius of the condensate. 



the condensate and vortices and strong bounds for the 
parameters. If dark matter in the galactic halo is com- 



posed of bosons in the form of a BEC with vortices then 
model independent arguments require a boson mass cen- 
tered around m ~ IQ-^o - IQ-^^ kg ~ IQ-^" - 10'^^ cV 
and a radius around 10^° — 10^^ m, which corresponds 
to roughly 1 — 100 vortices in the halo. Combining this 
with our solution, the scatting length is centered around 
a ~ 10~^° — 10"''^ m. Using our solution we also de- 
termined the size of the vortex core and calculated the 
critical rotational velocity for the nuclcation of vortices, 
showing that there exists large portions of parameter 
space where vortices would form. While these results are 
specific to the Andromeda galaxy, since the Andromeda 
galaxy is not atypical we expect them to hold more gen- 
erally. 
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